Critical wetting of a class of nonequilibrium interfaces: A computer simulation study 
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Critical wetting transitions under nonequilibrium conditions are studied numerically and ana- 
lytically by means of an interface-displacement model defined by a Kardar-Parisi-Zhang equation, 
plus some extra terms representing a limiting, short-ranged attractive wall. Its critical behavior is 
characterized in detail by providing a set of exponents for both the average height and the surface 
order-parameter in one dimension. The emerging picture is qualitatively and quantitatively different 
from recently reported mean-field predictions for the same problem. Evidence is shown that the 
presence of the attractive wall induces an anomalous scaling of the interface local slopes. 
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PACS numbers: 02.50.Ey,05.50.+q,64.60.-i 

I. INTRODUCTION 

Much scientific effort has gone into the study of equi- 
librium wetting since, in the late seventies, Cahn intro- 
duced the idea that it can be described as a phase tran- 
sition [1]. Among the various theoretical approaches de- 
veloped, interface displacement models have proved par- 
ticularly useful [2]. Within this perspective, the focus 
is on the interface that separates two coexisting (bulk) 
phases confined by a wall or substrate, and the wetting 
transition corresponds to the unbinding of the interface 
from the wall. This happens upon a rise of the tem- 
perature when the wall adsorbs preferentially one of the 
phases leading to a divergence of the thickness of the 
adsorbed layer. The dynamics of such an interface can 
be described at a coarse-grained level by the following 
continuum stochastic growth equation Q 



dh{x,t) 
dt 
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Here, /i(x, t) is the local height of the interface from the 
wall, the regions y > /i(x) and y < h{x) corresponding 
to the two bulk phases. D is the interfacial tension co- 
efficient and is a Gaussian white noise with zero mean 
and variance (t?(x, t)r]{x' ,t')) — 2a6{x ~ x')6{t — t') that 
mimics thermal fluctuations. V{h) accounts for the net 
interaction between the wall and the interface and its 
form depends on the nature of the forces between the 
particles in the bulk phases and with the wall, its rigor- 
ous derivation from microscopic Hamiltonians being far 
from trivial. If all the interactions are short-ranged, one 
may take in the limit of large h at phase coexistence 



V{h) 



dx 



(2) 



where T is the temperature Q . The amplitude c > is a 
repulsion whereas h{T) vanishes linearly with the mean- 



field wetting temperature, Tw^ as T — Tw^ and can rep- 
resent either an effective repulsion or attraction between 
the interface and the wall (see Fig. At sufficiently 
low temperatures, 6 < 0, the equilibrium thickness of the 
wetting layer as given by the stationary configurational 
average {h) is finite (pinned interface). This corresponds 
to an attractive potential (see Fig. [1]). As the temper- 
ature is raised the potential becomes less attractive and 
eventually, above a certain value 5 = 6^,, it no longer 
binds the interface and {h) diverges. Within mean-field 
approximation, ignoring spatial correlations, {h) follows 
from dV{h)/dh — 0, whereby one finds for an attractive 
wall (6 < 0) (h) — ln(— c/6) and, consequently, a critical 
wetting transition takes place as 6 — > bw = 0. Recently, 
effective short-ranged, equilibrium critical wetting show- 
ing mean-field like exponents seems to have been exper- 
imentally observed Q . 

Extensions of equilibrium, interface displacement mod- 
els to nonequilibrium conditions have only been recently 
addressed and constitute a topic of ongoing research ac- 
tivity [6j, 13] ■ Supplementing equation ([1]) with the most 
relevant nonequilibrium nonlinear term A(V/i)^ leads 
to a natural generalization of equation ([T]) that assumes 
that the velocity of the interface depends on its local- 
slope. 



dh{x, t) 
dt 



DV^h + X{\7hy + be-" + ce-^" + r]{x,t), (3) 



which is a Kardar-Parisi-Zhang (KPZ) interface Q in- 
teracting via a short-ranged potential with a wall. 
KPZ interfaces have a nonzero average velocity, v = 
A((Vft.)^), and hence steady-state interfaces move on av- 
erage thereby favoring one of the phases over the other. 
Wetting, on the other hand, by definition occurs at co- 
existence, i.e. zero average velocity of the free (no wall) 
interface. This agrees with the thermodynamic picture 
that at bulk coexistence any arbitrary fraction of the sys- 
tem may be in one phase, with the remainder in the other. 
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FIG. 1: (Color online) Mean-field binding potential for pos- 
itive and negative values of b. For 6 < the potential is 
attractive, exhibiting a well near the wall. 



Therefore, a constant Oc = —v needs to be included in 
equation ^ to study wetting transitions driven by the 
wall. The noncquilibrium analog of equilibrium critical 
wetting corresponds to the depinning transition at ac as 
b b~. Clearly, for A = and Ac = 0, the model reduces 
to the equilibrium one. In equilibrium, a constant force 
term, a, in the interfacial equation measures the devia- 
tion from bulk coexistence (i.e. it represents the chemical 
potential difference between the two phases) while here 
it plays a similar role by balancing the force exerted by 
the KPZ nonlinearity on the tilted parts of the interface, 
thereby guaranteeing that the average velocity of the free 
interface is zero. 

Owing to the lack of h ^ —h symmetry of the KPZ 
dynamics, it is necessary to specify either the relative 
position of the wall with respect to the interface, upper 
or lower, for a fixed sign of A, or reversely, the sign of 
A after a wall position has been arbitrarily chosen. In 
earlier studies of noncquilibrium (complete) wetting (see 
below), these two different physical situations lead to the 
existence of by-now-well-documented two different uni- 
versality classes, called multiplicative noise 1 (MNl) and 
multiplicative noise 2 (MN2) respectively [91]. Here, we 
take A < in Eq. ([3]) (i.e. the critical wetting coun- 
terpart of MNl) which has shown to have a much richer 
phenomenology than that of A > in other studies of 
nonequilibrium wetting . The analysis of the case 

A > will be tackled elsewhere. 

Thus defined, this model system, arguably the simplest 
nonequilibrium one, has served for the study of univer- 
sality issues in nonequilibrium wetting. In particular, by 
fixing 6 > i.e. in the presence of a repulsive wall, and 
letting a — > a~ (path 1 in figure [2]) nonequilibrium com- 
plete wetting transitions (MNl) were investigated. For 
h < hw (attractive wall) and varying a there is a rich 
phenomenology: the pinned and the depinned phases 
lose their stability at different values of a, giving rise 
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FIG. 2: (Color online) Schematic phase diagram for A < 
and a lower wall in the a — b plane (Eq. |4)|. The vertical 
line corresponds to the critical value a = ac, and the arrows 
denote the different types of transitions explained in the text. 
Path 1: complete wetting (upon approaching ac). Path 2: 
non-trivial depinning transition at a*{b). Path 3: First-order 
pinning transition. Path 4: Multicritical complete wetting. 
Path 5: Critical wetting. 



to a continuous depinning transition at a*(b) > a in the 
directed-percolation universality class (see Fig. [5J path 
2) [ll|, and a first-order phase-transition along path 3, 
at a = Qc- In the broad interval ac < a < a*{b) both 
phases coexist (see @ and references therein). Tricriti- 
cal behavior along path 4 was analyzed by Ginelli et al. 
p^ . and a preliminary study of nonequilibrium critical 
wetting (path 5) was presented in [isj . 

In this paper we investigate numerically and analyti- 
cally nonequilibrium critical wetting (path 5 in figure ^ 
as defined by equation ([3]). The focus will be on one- 
dimensional systems only. Higher system dimensionali- 
ties were studied in [To| by a mean-field analytic approx- 
imation to Eq. ^ which revealed the existence of three 
different regimes of scaling behavior. Their connection 
with our findings (or the lack of it) is discussed in the 
last section. 



II. MODELS AND OBSERVABLES 

Our continuous model is defined by the stochastic 
growth equation 



dt 



- \{Vhf + flc + be^'^ + ce"^'* + r^{x, t), 

(4) 

with, as explained above, Oc — — A((V/i)^). At the criti- 
cal wetting transition, i.e. as 6 approaches from below, 
the average stationary thickness of the wetting layer di- 
verges continuously as (h) ^ \b — 5^,1^'', where (3h is a 



critical exponent. At b = b^, {h{t)) ~ t^'^ for asymp- 
totically long times. Two other exponents we study are 
the dynamic and the correlation length exponents, z and 
V respectively, defined through their usual expressions 
^{t) ~ ~ |6 - bni]-", where ^ is the correlation 

length. They are related to the previous ones by the 
scaling form 9h — Pu/zv. 

Some of these exponents may be written in terms of 
known KPZ exponents. In particular, since s,tb — b^ the 
interface is asymptotically free, the dynamic exponent 
retains its one-dimensional free KPZ value z = 3/2 Q. 

Also, as we illustrate now. Oh is given by the expo- 
nent characterizing the growth of the interfacial width in 
the KPZ, W(t) - t^^^, and therefore. Oh ^ Ow = l/i 
in one dimension. This can be understood as follows: 
interfacial fluctuations are cutoff owing to the presence 
of the wall, as a result of which there is an effective 
fluctuation-induced repulsion between the wall and the 
interface. The latter can be estimated by noting that 
the wall makes itself felt when the mean interfacial sep- 
aration {h) is of the same order as the average extent 
of the interface fluctuations, 5h. From 5h ^ where 
C is the usual KPZ roughness exponent, and the defini- 
tion of V we find that the effective repulsion force has the 
form h~^^^^ and is therefore long-ranged in d = 1 where 
C = 1/2 > 0. Comparing this force with the determin- 
istic one in the Langevin equation ([?]), which is short- 
ranged, it is straightforward to conclude that in d = 1 
fluctuations dominate the unbinding of the interface. As 
fluctuations are governed, in the regime where the inter- 
face is asymptotically free, by the growth exponent of the 
KPZ, then Oh = 1/3. This result, has been verified in our 
computer simulations (see below). Note that this, as well 
as any other exponent computed exactly at the critical 
point, is path independent and, therefore, holds also for 
complete wetting. 

Of more interest is the behavior of the surface order- 
parameter defined as {e~^) or, equivalently as the density 
of local contacts between the interface and the wall. In- 
deed, considering equation ^ with D = —X, the change 
of variables h — — In n leads to 

^^^^^^^^ ^ DV^n-acn-bri^ -cn^ + nri{x,t), (5) 

which is a multiplicative noise Langevin equation for 
the surface order-parameter Q to be interpreted in the 
Stratonovich sense [l3]- For sufficiently low values of 
b < bw the interface remains pinned and the station- 
ary density of locally pinned sites at the wall is high 

^ !)• As the transition is approached (increasing b 
following path 5 in Fig. [2]), the stationary density of 
pinned segments goes to zero in a continuous manner as 
{n{b, t — oo)) ^ \b—bc\^'^ . Aib — b^,, the interface depins 
and therefore {h{t)) diverges and {n{t)) vanishes with the 
characteristic exponent {n{b — bw,t)) ^ t~^" . Obviously, 
as said before. On is common to paths 4 and 5 in figure [H 
and the value On = 0.5(1) has been reported previously 

. Here we focus here on the determination of the path 




FIG. 3: (Color online) Main: time decay of the surface order- 
parameter a.t q = q-w for simulations of the SSW model at 
system-sizes (from top to bottom) L = 2'^, 2*, 2^, 2^°, 2", 2^^, 
and 2^^" From the lowest curve S^^^ = 0.49(2). Inset: the 
crossover times to saturation as a function of the system-size 
lead to z^^^ = 1.4(1). The error bars indicate the estimated 
uncertainties in the crossing point of fits of the initial decay 
and the saturating behavior to straight lines, on the log-log 
plot. 



dependent exponents f3h, Pn, and v. Before proceeding 
further, we refer the reader to [lB| for a detailed anal- 
ysis of the critical behavior of the above observables in 
equilibrium critical wetting. 

To study numerically the model defined above, ow- 
ing to well documented numerical instability problems 
|la |. it is more convenient to integrate numerically Eq. 
^ than the equivalent form Eq. (j4|). Equation ^ 
can be efficiently integrated by means of a recently in- 
troduced split-step scheme specifically designed to deal 
with Langevin equations with non-additive noise p^j . 
Setting = — A = 0.1,(T = lwe can use the result 
flc = 0.143668(3) obtained from previous investigations 
of the analogous nonequilibrium complete-wetting tran- 
sition [l8|. As will be illustrated below, estimates of the 
critical exponents are severely hindered by uncertainties 
in the value of Uc- 

To circumvent this problem and to confirm universality 
we have carried out simulations of a discrete interfacial 
growth model which (i) in the absence of walls is known 
to belong to the KPZ universality class, (ii) has been suc- 
cessfully used in nonequilibrium complete wetting anal- 
yses and, most importantly, (iii) allows for an exact 
determination of the velocity of the free interface, and 
therefore permits to extract the critical exponents with 
good accuracy. To be more specific, we consider a single 
step plus wall model (SSW) defined as follows. At time t 
interface positions above sites i of a one-dimensional line 
of length L are given by integer height variables ht{i), sat- 
isfying the solid- on- solid constrain \ht{i) — ht{i + \ )\ — 1. 
Initially, we take ho{2i) — and /io(2i -I- 1) = 1. New 
height configurations are generated by choosing at ran- 
dom a site i and growing it to ht{i) — > ht{i) -I- 2 if and 



only if a local minimum existed at i. It can be shown that 
this rule generates a KPZ-like interface with A = — 1/2 
moving with an asyrnptotic long-time, average velocity 
V = {1 + L)/{2L) [20]- Aditionally, taking advantage of 
the exact knowledge of v, the interface is globally pulled 
down by one unit every L/v growth trials, in such a way 
that the interface has zero average velocity. A wall at 
/i = is then introduced by precluding the interface from 
overtaking the wall that is behind it. This is achieved by 
implementing the previous global, downwards movement 
as ht{i) = \ht{i) — l\ for each i. Finally, in order to imple- 
ment an attractive wall the growing rates at the bottom 
layer are reduced from 1 to 1 — g with < q < 1. The pa- 
rameter q represents the short-ranged attraction exerted 
by the wall on the interface. If g = 1 no growth is possi- 
ble at the local minima located at the wall and hence an 
interface at the wall does not move, whereas if q = the 
short-range attraction is switched off. 

In all, the SSW algorithm is as follows: (i) A site i 
is randomly chosen and grown from h{i) to h(i) + 2 if a 
local minimum exists at i {h{i + 1) -|- h(i — 1) — 2h(i) ~ 
2). This is done with probability 1 if h{i) > 0, or with 
probability 1 — q if h{i) = 0, with < q < 1. No action 
is taken if i does not correspond to a minimum. Time 
is increased by l/v after L of such attempts, (ii) Every 
L/v growth trials h{i) = \h{i) — 1| for each i. Given that 
L/v is generally not an integer, this is done by using 
[L/v\ with probability L/v — [L/v\ and [L/v\ +1 with 
probability [L/vJ -1-1 — L/u, where [-J denotes the integer 
part. Periodic boundary conditions are imposed. By 
tuning the parameter q a critical wetting transition is 
observed. 



III. NUMERICAL RESULTS 

We next summarize our main findings for both the dis- 
crete model (SSW) and the stochastic differential equa- 
tion dH) (SDE). 

To determine the critical points, q^ (for the SSW) 
and (for the SDE), we take a system-size as large 
as possible {L = 2^^ here), plot the order-parameter 
(n(t)) versus i in a double logarithmic scale, and look 
for the separatrix between curves converging to a con- 
stant value and those bending downward (not shown). 
In this way, the critical values g„ = 0.4445(5) (SSW) 
and bw = —1.0815(9) (SDE) are determined. In figure 
13] results are shown in log-log for the time decay of the 
order-parameter {n{t)) at the critical point of the SSW 
model for different system-sizes ranging from L = 2^" to 
L = 2^^. From the slope of a straight-line fit to the lowest 
curve one finds ^^sw ^ 0.49(2) {e^^^ = 0.50(5) for the 
SDE; not shown), where the error is computed by com- 
paring the slopes corresponding to the upper and lower 
bounds of q^ (bw)- The inset of Fig. [3] shows estimates 
of the crossover times, tx (L), from time decay to satura- 
tion as a function of the system-size. Using tx{L) 
we find the temporal exponent values z^^w _ 1.4(1), 




FIG. 4: (Color online) Main: The SSW saturation values of 
(n) at criticality for several system-sizes provide the exponent 
ratio Pn/i^ ~ 0.74(1). Inset: the same analysis for the SDE 
yields 0.6(1). Note the difference in the magnitude of error- 
bars, which correspond to three standard deviations of the 
mean saturation values. 



^.SDE _ 1.3(2) (not shown), which are in agreement with 
those previously reported in [l3| following path 4 of Fig. 
[21 and compatible with the theoretical considerations de- 
scribed above. 

From the scaling of the saturation value at the critical 
point, {usatiq = q-w)) for different system-sizes one can 
determine (see figure S]) /3„/z/ = 0.74(1) {(3n/y = 0.6(1) 
for the SDE). A direct estimation of /3„ is also possible by 
measuring the order-parameter stationary-value for the 
largest available system-size upon approaching the criti- 
cal point from below. This is shown in figure [5] again for 
both the SSW and the SDE. We find /3,f ^ ^ 1.50(9) and 
(3^^ — 1.46(6) which, along with the obtained values for 
I3n/v, yields v^^^ = 2.0(2) and u^^^ = 2.4(5) (note the 
relatively large errorbar in this latter case). These values 
supersede the early estimate /9„ = 1.2 given in [l^ 

The scaling properties of the mean interfacial sepa- 
ration {h) can be determined analogously (we show re- 
sults only for the SSW model). The time growth of 
{h{t)) for the largest system available L = 2P yields 
^h^^ = 0.35(2) (see Fig. [S]), in reasonable agreement 
with the expected value Qh = 1/3. The average satu- 
ration values are plotted in log-log as a function of the 
system-size in the inset of figure ^ . From them we esti- 
mate (ih/v — 0.52(4) (compatible with KPZ scaling, as 
C — Ph/v = 1/2). Estimations of z can be analogously 
obtained, but they are rather noisy because of the uncer- 
tainty in determining the saturation value. Also, a direct 
estimation of Ph can be obtained for the largest available 
size. This leads to /3f ^ = 0.9(1) and = 1.0(1) 

(see figure [7]). Using the values of Ph/i^ and (3^ a third 
estimate for ly — 1.8(3) is obtained. 

Finally, we have confirmed that for both models at 
the tricritical point the width, W, of the interface grows 
with time with an exponent compatible with that of the 




FIG. 5: (Color online) Main: log- log plot of the satura- 
tion values of (n) for L = 2^^ vs the distance to the critical 
point gives a direct estimation of (5^^ = 1.50(9) and (inset) 
PI^^ = 1.46(6). Error bars as in Fig. H 




FIG. 7: (Color online) Main: the scaling of the saturation 
value of (hsat) yields 13^^^ — 0.9(1) for data collected from 
the SSW model. Inset: A value (3^°^ = 1.0(1) results for data 
collected from the SDE. Error bars as in Fig. |4] 




FIG. 6: (Color online) Main: log-log plot of the average 
distance to the wall vs time, yielding 0^^^ = 0.35(2). System- 
sizes are (from top to bottom) L = 2^'^, 2^^, 2", 2^^^, 2^^, 2®, 
and 2^. Inset: finite-size scaling of {hsatil = q-w)) indicating 
Ph/i' = 0.52(4). AU values are for the SSW model. Error 
bars as in Fig. 3] 



KPZ, W{t) with Ow = 1/3, and saturates in 

finite system-sizes to a non- vanishing value given by W ^ 
L-f^^l^ ~ L-i/2 as in the KPZ. 

In summary, the exponent values determined above ap- 
pear to be compatible with the set of simple rational 
numbers 6'„ = 1/2, z = 3/2, /3„ = 3/2, i/ = 2,6*/, = 1/3, 
and /3/i — 1. For the sake of comparison, the critical ex- 
ponents at equilibrium critical wetting are On = 1/4, z = 
2,/3„ = 1,!/ = 2,61,, = 1/4 and /3ft = 1 [2l|. The coinci- 
dence in the numerical values of two of them (the path 
dependent ones, (5h, and v) is somehow intriguing. 



IV. THEORETICAL CONSIDERATIONS 

It is instructive to compare the above results with 
those obtained at high system dimensionalities from a 
self-consistent, mean- field approximation to equation 
p^ . In that approximation, a sequence of three scal- 
ing regimes was reported to exist depending on the rela- 
tive importance of the noise strength as compared to the 
spatial coupling. The first two have a Gaussian charac- 
ter, while in the third regime, i.e. the strong noise one 
(in which the wetting temperature is shifted away from 
zero) all moments, = (n''), for k > 2D/cr'^ + 1 scale 
with the same exponent, while simple scaling 
is obtained for k < 2D/a^ + 1. This is a rather cu- 
rious type of anomalous scaling not very different from 
that observed in analogous mean-field approximations for 
non-equilibrium complete wetting 22]. However, we have 
verified numerically that in the one-dimensional system 
moments of arbitrary order scale as (n) itself, as happens 
in other Langevin equations with multiplicative noise 0] . 
This is a consequence of the large fluctuations occurring 
at d = 1, notwithstanding which it is possible to make 
some analytic predictions, as we now discuss. 

First, we explicitly show how the critical value of b is 
depressed from its mean-field value 6^, = to 6^ < once 
fluctuations are included. Taking spatial and noise aver- 
ages in the stationary state of equation (|3|) for a generic 
value of o, and denoting the average squared slope of the 
interface by = ((V/i)^), we obtain 

a-s^ + b{n) + (n^) = (6) 

where, without loss of generality we have taken D — 
A = — 1 and c = 1. At coexistence Uc — s^, and hence 
— + b{n) + {m?) = 0. Additionally, in the pinned 
phase (i.e. (n) ^ 0) < and therefore equation ([6]) 
has a solution only if 6 < 0. Still, 6^, could be 0, but 
by noting that on approaching the tricritical point along 
path 5 the moments of (n) scale in the same way, {n) ~ 
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{n?) ~ \h — hw\^^ , it is easy to see that 6„ 9^ because 
otherwise the negative term b{n) would be subdominant 
in comparison with the positive (ri^), and no solution 
could exist for small h. Clearly, by continuity hw cannot 
be positive and hence 6^ < ensues. 

It is illuminating to show how a similar reasoning leads 
to a more predictive analysis when applied to the anal- 
ogous complete wetting transition (path 1 in Fig. ^ 
known to be in the MNl class. In that case, the con- 
dition a — ac -\- — + h{n) =0 must be satisfied, with 
the negative term 5a = a — balancing the two posi- 
tive terms sj. — and h{n) ~ jJaj''" in the pinned phase. 
This implies that /3„ > 1 if a solution for (n) is to exist 
for small 5a Furthermore, by noticing that if KPZ 
scaling is applicable, then — ~ ^2(i-c) ^23*], where C 
is the roughness exponent of a free KPZ. Recalling that 
C = 1/2 in Id, immediately entails v = \ for complete 
wetting [7]. 

On the contrary, we have not been able to derive the 
value of V for nonequilibrium critical wetting, nor does 
it seem immediate that v = \ along path 4 as found 
numerically in [T2|. The main difference with the com- 
plete wetting case appears to rest on the behavior of the 
slopes . According to our measurements for equation 
(HI), si- r^\b- for critical wetting (path 5) 

and si — ^ \a — ad^''^*'-^-' for tricritical complete wet- 
ting (path 4), pointing to an anomalous scaling of the 
slopes, i.e. = ((V/i)^) does not have the same scal- 
ing dimension as as happens in the complete 
wetting case. 

That si — scales along path 4 with an exponent less 
than unity could have been anticipated from the condi- 
tion \5a\ — bw\5a\^'^ — s^. — s^ at the tricritical point, 
implying s^ — s"^ ^ \5a\^'' with /3s < 1- 

The exponent value /3„ = 0.74(5) reported in [l2| 
for the multicritical complete wetting transition along 
path 4 (confirmed in our own measurements), along with 
/3„ w 1.5 for critical wetting as referred above, indicates 
that the anomalous scaling of the slopes is ultimately 
controlled by /3„ in both cases, rather than by the rough- 
ness exponent of the KPZ as in the complete wetting 
case. Notice the constancy along path 4 and 5 of the 
ratio (3n/v ~ 0.75, which is a property of the tricritical 
point and therefore must be independent of the path. 

The fact that the slopes acquire a scaling not directly 
derivable from the free KPZ can be argued to be a con- 
sequence of the effect of the potential well on the inter- 
face. Figure [5] shows snapshots of configurations h{x) 
that result from solving the SDE in the stationary state 
for different parameters that correspond to approaching 
the critical points from different paths: 

• Panel A for (6 = -1.2 < bw,a = -0.25 < Cc); 
corresponds a situation slightly below the complete 
wetting (MNl) transition (path 1 of Fig. 

• Panel B for (6 — h^^a = —0.153 < Ac), slightly 
below the multicritical complete wetting transition 
(path 4). 



• Panel C for (& = — 1.3 < 6^,, a = Cc), slightly below 
the critical wetting transition (path 5). 

The distances to the transition points are chosen so 
that (h) « 5 in all cases. Note the clear qualitative 
difference between panel A, for which standard scaling 
holds, and the rest. Observe that in panel C, the inter- 
face consists of patches of essentially free KPZ interfaces 
separated by regions of sites pinned by the potential well. 
It then seems plausible to conclude that, following path 5, 
regions locally trapped within the potential well develop 
slopes different from that of a free KPZ (indeed, rough- 
ness is severely restricted within the potential well) even 
at points arbitrarily close to the tricritical point, thereby 
inducing an anomalous scaling controlled by /?„ . In panel 
B, i.e. upon approaching the tricritical point along path 
4, the effect of the bounding wall is less apparent as the 
potential well is marginally disappearing at b ~ h^, but 
a similar effect, induced by the potential shape, should 
be at work in this borderline case. 



V. SUMMARY AND DISCUSSION 

We have investigated the universal properties of 
nonequilibrium critical wetting transitions in one spatial 
dimension. For that, we study effective interfacial mod- 
els in the KPZ universality class with A < bounded by 
a lower wall, and determine the average height {h) and 
the surface order-parameter (n) = (e"''). 

A scaling analysis leads to the prediction z = 3/2 and 
a time behavior of the average height governed by the 
growth exponent of the free KPZ, i.e. {h{t)) ~ t^/^. 
These results have been verified numerically. Other expo- 
nents have been computed from extensive numerical sim- 
ulations of the Langevin equation and a discrete model 
in the same universality class that enables a more precise 
numerical analysis, as a result of which we find ~ 2, 
/3/i ~ 1, and /3„ ~ 3/2, suggesting that actually the ex- 
ponents take rational values (see table I). Interestingly 
enough, the first two values agree with those of equilib- 
rium critical wetting. 

Simple analytical arguments allow us to show that the 
critical value of the control parameter b is depressed by 
fluctuations from its mean-field value fe^ = to 6^, < 0. 
We have also shown that in critical wetting, as well as 
in multicritical complete wetting, the average interface 
slopes do show anomalous scaling, not controlled by the 
free KPZ equation: local regions pinned by the binding 
potential generate anomalous scaling. 

We have not been able to predict the value of v us- 
ing the same analytical considerations that yield v — \ 
in the complete wetting (MNl) case, nor is it trivial to 
obtain = 1 for multicritical complete wetting. It is nev- 
ertheless possible to conclude that in this latter case the 
exponent governing the scaling of the interface slopes is 
less than unity, — ~ |a — ad , in agreement with the 
value (3s ~ 0.75 obtained. There is numeric evidence that 
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FIG. 8: (Color online) Snapshots of configurations h{x) as 
results from solving the SDE in the stationary state for (A) 
b = —1.2 < bu,,a = —0.25 < Oc (complete wetting), (B) 
b = bw,a = —0.153 < Uc (multicritical complete wetting), 
and (C) b — —1.3, a = ac (critical wetting). The distances 
to the transition points are chosen so that (h) « 5 in all 
cases. In panels (B) and (C) patches of depinned interfaces 
are observed where n ~ 0. 



TABLE I: Summary of the critical exponents for nonequi- 
librium, critical wetting transitions with short-range forces. 
Results are shown for the discrete model SSW, the stochas- 
tic differential equation SDE, and from other sources when 
available. 



Exponent 


SSW 


SDE 


Others 


z 


1.4(1) 


1.3(2) 


1.5(1) [11] 


On 


0.49(2) 


0.50(5) 


0.50(1) [11] 


9h 


0.35(2) 


0.33(1) 




V 


2.0(2) 


2.4(5) 






1.50(9) 


1.46(6) 




I3h 


0.9(1) 


1.0(1) 





the slopes actually scale with Ps — /3„ for both critical 
wetting and multicritical complete wetting (paths 5 and 
4 of Fig. [21 respectively), a result that violates naive scal- 
ing. In effect, dimensional analysis demands that — 
scales as \a — ac\^~''\ C being the roughness exponent 
of the free KPZ [23[ . This is obeyed at the complete wet- 
ting transition, and indeed was used to derive j/ = 1 in 
MNl Q, but does not hold for critical wetting nor for 
multicritical complete wetting. 

The cause of the deviation from standard scaling can 



be sought in the interfacial profiles shown in figure[8] For 
critical and multicritical wetting, patches of free KPZ in- 
terfaces (n(a;) « 0) are separated by regions of sites that 
lie in the potential well, plausibly hindering the standard 
KPZ scaling from setting in gradually as the tricritical 
point is approached. This is in contrast to typical interfa- 
cial profiles for noncquilibrium complete wetting (MNl) 
where no potential well is at play. 

According to a recent self-consistent, mean-field ap- 
proximation to equation ([5]) three different scaling 
regimes of critical behavior for the surface order- 
parameter can be distinguished The first two are of 
Gaussian type, while the third one is a highly nontrivial 
strong-fluctuating regime. This rich structure is com- 
pletely washed out by fluctuations in one-dimensional 
systems, where a unique and universal scaling regime 
emerges. Moreover, we have verified numerically that in 
the one-dimensional system moments of arbitrary order 
scale as (n) itself, as happens in other Langevin equations 
with multiplicative noise 0. The fact that the numer- 
ical values of the exponents are changed is not surpris- 
ing at all, given the presence of severe fluctuations in 
one-dimension but what is more striking is that out of 
the three regimes appearing in the mean-field approach, 
only one survives. A possible explanation for such an 
abrupt change might come from a recent claim that the 
strong-coupling renormalization group fixed point of the 
KPZ dyriamics is essentially different above and below 
d = 2 [2J]. This point remains to be further studied, 
as well as some other aspects of noncquilibrium critical 
wetting including its subtle relation to its equilibrium 
counterpart (signaled by the coincidence of some expo- 
nents), and the possible existence of various nonuniversal 
scaling regimes in higher dimensions. It also stands as a 
main experimental challenge to observe in the laboratory 
the phenomenology reported on here and in previous the- 
oretical works of noncquilibrium wetting. 
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